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ABSTRACT 

In  this  paper  we  provide  an  accurate  analysis  of  the  performance  of  a  random-carrier  (RC) 
code-division  multiple-access  (CDMA)  scheme  recently  introduced  for  use  in  high-capacity  optical 
networks.  According  to  this  scheme  coherent  optical  techniques  are  employed  to  exploit  the  huge 
bandwidth  of  single-mode  optical  fibers  and  are  coupled  with  spread-spectrum  direct-sequence 
modulation  in  order  to  mitigate  the  interference  from  other  signals  due  to  the  frequency  overlap 
caused  by  the  instability  of  the  carrier  frequency  of  the  laser,  or  to  the  mistakes  in  the  frequency 
coordination  and  assignment. 

The  average  bit  error  probability  of  this  multiple-access  scheme  is  evaluated  by  using  the  char¬ 
acteristic  function  of  the  other-user  interference  at  the  output  of  the  matched  optical  filter.  Both 
phase  noise  and  thermal  noise  are  taken  into  account  in  the  computation.  Time-synchronous  as 
well  as  asynchronous  systems  are  analyzed  in  this  context.  Binary  phase-shift-keying  (BPSK)  and 
on-off-keying  (OOK)  data  modulation  schemes  are  considered.  The  analysis  is  valid  for  arbitrary 
values  of  the  spreading  gain  and  the  number  of  interfering  users.  The  performance  evaluation  of  RC 
CDMA  establishes  the  potential  advantage  in  employing  hybrids  of  wavelength-division  multiple- 
access  (WDMA)  and  CDMA  to  combat  inter-carrier  interference  in  dense  WDMA  systems. 


This  research  was  supported  in  part  by  the  Office  of  Naval  Research  under  contract  N00014- 
89- J- 1375  and  in  part  by  the  Systems  Research  Center  at  the  University  of  Maryland,  College 
Park,  through  the  National  Science  Foundation’s  Engineering  Research  Centers  Program,  NSF 
CDR  8803012. 


1.  Introduction 

The  random- carrier  (RC)  code-division  multiple- access  (CDMA)  scheme  was  first  in¬ 
troduced  in  [l]  for  use  in  high-capacity  optical  networks.  According  to  this  scheme,  coher¬ 
ent  optical  techniques  are  employed  to  exploit  the  huge  bandwidth  (tens  of  thousands  of 
GHz)  of  single-mode  optical  fibers.  The  inherent  instabilities  of  present-day  semiconductor 
lasers  are  circumvented  by  coupling  the  optical  multiple-access  system,  which  is  assumed 
to  place  randomly  the  modulated  carriers  in  the  available  optical  band,  with  CDMA.  In 
particular,  spread-spectrum  direct-sequence  modulation  is  employed  in  order  to  mitigate 
the  interference  from  other  signals,  that  is  due  to  the  frequency  overlap  caused  by  the 
instability  of  the  carrier  frequency  of  the  laser. 

Analysis  of  RC  CDMA  provides  the  basis  for  exploiting  the  potential  advantage  of 
employing  hybrids  of  wavelength-division  multiple-access  (WDMA)  and  CDMA  to  combat 
inter-carrier  interference  in  dense  WDMA  schemes.  That  hybrid  scheme  used  WDMA  for 
providing  multiple-access  capability  and  CDMA  for  protection  against  laser-frequency  in¬ 
stabilities  and  mistakes  in  the  frequency  coordination  and  assignment.  Such  a  combination 
of  the  best  features  of  both  WDMA  and  CDMA  wa.s  recently  proposed  in  [2]. 

We  complete,  extend,  and  validate  the  work  of  [1]  by  providing  a  more  accurate  evalu¬ 
ation  of  the  performance  of  the  RC  CDMA  scheme,  without  making  all  the  approximations 
and  limiting  assumptions  made  there.  Here  are  the  specific  contributions  of  our  work:  (a) 
The  performance  measure  considered  is  the  average  probability  of  a  bit  error  at  the  out¬ 
put  of  the  RC  CDMA  system’s  optical  matched  filter  and  is  evaluated  via  an  approach 
that  can  achieve  any  desirable  accuracy,  (b)  Our  analysis  is  valid  for  an  arbitrary  spread¬ 
ing  gain  (number  of  chips  per  bit)  and  an  arbitrary  number  of  interfering  users,  (c)  In 
this  evaluation,  we  also  take  into  account  the  effects  of  phase  noise  of  the  lasers  and  shot 
noise;  the  former  is  modeled  as  a  Brownian  motion  process  and  the  latter  as  additive  white 
Gaussian  noise  (AWGN).  (d)  Not  only  chip-synchronous,  but  also  asynchronous  multi-user 
systems  are  analyzed;  thus  the  effect  of  time  delays  between  the  various  users  are  taken 
into  account,  (e)  Besides  PSK,  on-off-keying  (00K)  data  modulation  is  also  considered. 

For  the  evaluation  of  the  bit  error  px'obability,  we  use  the  characteristic-function 
method  introduced  in  [3]  for  radio- frequency  (RF)  CDMA  communication  systems.  The 
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evaluation  of  error  probability  in  [3]  was  carried  out  for  arbitrary  deterministic  signature 
sequences.  In  [4],  the  results  of  [3]  were  extended  to  CDMA  systems  employing  random 
signature  sequences  (i.e.,  i.i.d.  sequences  that  assume  the  values  +1  and  —1  with  equal 
probability  and  are  mutually  independent  for  different  users).  In  this  work,  we  analyze 
optical  RC  CDMA  systems  with  random  signature  sequences.  In  our  analysis,  BPSIv  mod¬ 
ulation,  as  well  as  00K  modulation  is  used  to  modulate  the  data  bit  stream,  while  M-ary 
PSI\  modulation  is  employed  for  the  signature  sequence  stream.  Electro-optical  phase 
modulators  [6]  make  these  features  feasible. 

We  evaluate  the  average  bit  error  probability  by  averaging  over  the  data  streams, 
signature  sequences,  carrier  frequencies,  phases  of  the  interfering  users,  and  time  delays 
(for  asynchronous  systems).  This  is  accomplished  by  computing  the  characteristic  function 
of  the  interference  due  to  the  other  users  at  the  output  of  the  optical  filter  that  is  matched 
to  a  particular  signal.  The  accuracy  of  this  computational  technique  can  be  completely 
controlled  by  the  user  and  is  determined  by  the  accuracy  of  the  integration  routines  invoked: 
as  the  the  number  of  points  in  the  integration  rule  increases,  the  required  computer  CPU 
time  increases.  Any  desirable  accuracy  can  be  attained  with  the  help  of  this  technique. 

In  the  calculation  of  the  characteristic  function  of  the  multiuser  interference  the  Cen¬ 
tral  Limit  Theorem  (CLT)  is  used  once;  however,  its  use  is  properly  justified  and  the  results 
are  compared  to  those  obtained  via  another  approach  that  makes  no  use  of  the  CLT. 

In  our  analysis,  the  characteristic  function  method  is  actually  applied  to  more  general 
situations  than  those  considered  in  [3] -[5],  which  deal  with  other-user  interference  and 
multipath  interference  and  have  a  symmetric  probability  density  function  (pdf).  In  our 
work,  the  characteristic  function  technique  used  for  the  evaluation  of  the  average  error 
probability  of  systems  with  additive  interference  is  extended  to  the  case  of  interference 
with  non-symmetric  pdfs  (the  case  when  the  interference  does  not  assume  positive  and 
negative  values  with  equal  probability). 

Another  useful  product  of  our  analysis  of  RC  CDMA  systems  is  an  exact  expression 
for  the  error  probability  of  a  single-user  coherent  optical  system  disturbed  by  phase  noise 
and  AWGN. 

We  provide  numerical  results  describing  the  performance  of  the  RC  CDMA  systems. 
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The  tradeoffs  between  the  various  system  parameters  are  illustrated  and  interpreted.  Be¬ 
sides  the  average  error  probability  as  a  function  of  system  parameters,  the  maximum 
number  of  users  that  can  be  supported  with  this  scheme  at  a  given  error  probability  is  also 
obtained.  Moreover,  the  performance  of  a  single  user  coherent  optical  system  disturbed 
by  phase  noise  and  AWGN  is  analyzed  in  detail. 

Our  analysis  of  the  RC  CDMA  system  paves  the  way  for  accurately  evaluating  the 
performance  of  dense  WDMA  systems  and  of  hybrid  WDMA/CDMA  systems,  and  for 
comparing  the  performance  of  such  systems  in  terms  of  average  error  probability  and 
multiple-access  capability  (i.e.,  the  maximum  number  of  simultaneous  users  that  can  be 
supported  at  a  prespecified  error  probability)  for  realistic  environments  characterized  by 
asynchronous  users,  laser  phase  noise,  shot  (AWGN)  noise,  and  frequency  instability  of 
lasers.  In  Chapter  3,  we  provide  the  analysis  of  pure  WDMA  and  hybrid  WDMA/CDMA 
systems  and  a  comparison  of  their  performance. 

This  paper  is  organized  as  follows.  In  Section  2,  the  model  and  the  receiver  structure 
are  described.  In  Section  3,  the  single  user  performance  for  the  BPSK  and  OOIv  modulation 
is  obtained.  Section  4  extends  the  analysis  to  the  multiuser  case,  in  which  K  active  users 
share  a  common  optical  channel.  The  average  bit  error  probability  for  the  intended  user 
is  obtained  by  using  the  characteristic  functions  of  the  interference  and  AWGN.  Section 
5  contains  computation  of  the  characteristic  function  of  other-user  interference.  Based  on 
two  different  sets  of  assumptions,  this  function  is  obtained  for  both  BPSK  and  OOIv  mod¬ 
ulation.  In  Section  6,  the  pdf  of  a  useful  random  variable,  which  plays  an  important  role 
in  the  analysis  of  both  single  user  and  multiuser  systems,  is  estimated.  Section  7  extends 
the  analysis  to  the  asynchronous  system,  in  which  random  delays  are  introduced  to  the 
signals  of  the  interfering  users.  Section  8  presents  the  numerical  results  and  comparisons 
of  different  systems.  Finally,  Section  9  contains  the  summary. 

2.  Model 

K  high  data  rate  users  share  a  common  optical  channel  in  a  multi-access  fashion. 
These  users  are  scattered  in  an  optical  bandwidth  as  big  as  W  =  10 THz.  Due  to  the 
frequency  instabilities  of  lasers,  each  carrier  frequency  wanders  around  its  designated  fre- 
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quency.  We  assume  that  each  carrier  frequency  is  randomly  distributed  in  the  total  optical 
bandwidth  W .  The  transmitted  optical  signal  is  denoted  by  5(f),  which  is  a  complex  signal 
as 

I< 

S(t )  =  5Z  VP  bm(t  -  Tm)am(t  -  Tm)e^t-rm)+9mw]  (2.3.) 

m— 1 

where  associated  parameters  are  as  follows: 

•  P  is  the  transmitted  signal  power  of  each  user 

•  brn(t)  is  the  data  stream  of  the  m-th  user  given  by 

OO 

bm(t)=  Y  b[m)p(t-nT) 

n  —  —  oo 

where  b denotes  the  n-th  bit  of  the  m-th  user;  — 1, 1}  for  BPSIv  modulation 

and  bn  V{0, 1}  for  OOIv  modulation:  pit)  is  a  pulse  of  unit  amplitude  in  [0,  T]. 

•  om(f)  is  the  addressing  function  or  signature  sequence  stream  used  by  user  m.  That 
is 

OO 

am(t)  =  e^f)  =  Y  e'^mn h(t  —  nTc) 

n~  — oo 

where  h(t)  is  a  pulse  of  unit  amplitude  in  [0,  Tc]  and  Tc  =  is  the  chip  duration, 
where  N  is  the  number  of  chips  per  bit.  is  a  phase  taking  values  in  [0,  2tt] . 

•  u>m  is  the  carrier  on  which  the  m-th  signal  is  sent.  This  value  is  randomly  chosen  by 
the  transmitter  laser  for  (RC)  CDMA,  or  is  preassigned  for  hybrid  WDMA/CDMA. 

•  Qm(t)  is  the  phase  noise  associated  with  the  mth  transmitter  laser,  which  is  a  Brownian 
motion  process  with  Lorentzian  bandwidth  /3.  The  mean  of  this  process  is  zero  and 
the  variance  is  2ir (3t. 

•  Tm  is  the  m-th  time  delay  which  is  a  uniform  random  variable  distributed  in  [0,  T]. 
For  the  synchronous  system,  this  time  delay  is  zero. 

At  the  &-th  receiver,  the  optical  signal  5(f)  is  first  despread  by  a*k(t),  which  is  the 
complex  conjugate  of  a.k(t )  ,  and  then  homodyne-detected  for  the  transmitted  signal  from 
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user  k  (see  Fig.  1).  The  output  of  the  photodetector  is 
r(t)  =  VP  bk{t)elA9k{t)  +n{t)+ 

t 

Vp  _  Tm^ei[u'rn(t-Tm)+<t>rn(t—Trn)-<j>k(t-Trn)+A9rn{t)} 

m 


0  0\ 


E =  E 

m  m  =  l 

m  ^  k 

I  A 

^772  —  ^  A: 

&0k{t)  =  6k(t)  -  $L(t) 

Aem(t)  =  em(t)-eL(t) 

where  #£,(£)  is  the  phase  noise  of  the  local  laser  and  n{t)  is  the  complex  AWGN  process 
with  double-sided  spectral  density  -qr  .  The  receiver  used  is  a  correlation  receiver,  which 
is  optimum  for  the  single  user  case  with  no  phase  noise  (see  Fig.  2).  The  performance  of 
this  suboptimum  receiver  in  the  presense  of  phase  noise  and  AWGN  is  obtained  for  the 
single  user  and  multiuser  situations  in  the  subsequent  sections. 
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3.  Single  User  Analysis 

Here  the  performance  of  the  system  described  in  Section  2,  is  evaluated  in  terms  of 
BER  for  the  single-user  case.  The  real  part  of  the  output  of  the  integrator  is  denoted  by 
Y,  which  is 

Y  =  b(Qk)xVP  +  i}s/P  (3.1) 


where 


(3.2) 


and  77  is  a  zero  mean  Gaussian  random  variable  of  variance  No/'2PT.  According  to  Ap¬ 
pendix  A,  the  probability  of  error  Pe  for  BPSK  modulation  is  obtained  as 


^(/>  +  X ) 


+  -$  ( p-X 


(3.3) 


where  Q(-)  and  $(•)  are  related  to  the  standard  normal  distribution  as  follows: 


The  overlines  in  (3.3)  indicate  expectation  with  respect  to  the  random  variable  X.  Simi¬ 
larly,  for  00K  modulation.  Pe  is  obtained  as 


+  I4 


Up  -  X) 


(3.4) 


It  is  easy  to  show  that  the  optimum  value  of  p  for  (3.3)  and  (3.4)  is  0  and  1/2.  respectively. 
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4.  Multiuser  Analysis 

Here  the  performance  of  the  system  of  Section  2.1  is  evaluated  in  terms  of  BER  for 
the  multiuser  case.  The  synchronous  case  (rm  =  0)  is  considered  first;  the  asynchronous 
case  is  studied  in  a  separate  section.  The  output  of  the  integrator  V  is 

V  =  b{0k)VPx  +  \/pJ2im  (4.1) 

m 


where 


and 


T 

i  /  e‘AI>Mdt 

T  Jo 


-  b(m)  ■  — 

—  Un  — 


dt 


Ts/p 


n(t)clt. 


(4.2) 

(4.3) 


(4.4) 


In  order  to  calculate  the  integral  in  (4.3)  we  make  some  assumptions.  First,  we  assume 
that  the  phase  addressing  functions  4>m(t)  and  <j>k(t)  in  the  chip  interval  ({n  —  1)  Tr ,  nTc ] 
take  the  values  <J)mn  and  <f>kn  in  [0,27r],  respectively.  Second,  for  large  values  of  N.  we 
assume  that 


A 6m(t)  =  A 9m(nTc)  =  dmn,  for  ( n  -  1  )TC  <  t  <  nTc  (4.5) 

where  is  a  zero-mean  Gaussian  random  variable  of  variance  4v r'3nTc  .  Under  these 
assumptions, 


2  m  — 


,(m) 

J_0 _ 

N 


-sine 


u'mTc 


N 

E 

n=  1 


ij^mn  —  ©fcn+#  m  n  Jru>'m  (  n~  1/2)  rc] 


(4.6) 


By  taking  the  real  part  of  V  in  (4.1),  we  have  ' 


Y  =  b[k)  -X  +  Y,  tm  +  V 


m 


(4.7) 


where  X  and  77  are  defined  in  (3.1)  and  (3.2)  and  Im  is 


'XT, 


N 


-sine 


N 


n—  1 


cos  (A', 


(4.S) 


where 

Xmn  =<  (f>mn  —  <j>kn  +  #mn  +  Xm  ( 77  —  1/2 )TC  >  .  (4.9) 

In  (4.9),  <  •  >  represents  [•]  mod.2iv  .  To  evaluate  the  bit  error  probability  let  us  rewrite 
(4.7)  as 

Y  =  ( b{0k)  ■  X  +  I  +  r/)  >/P  (4.10) 

where 


(4.11) 


According  to  Appendix  A,  the  average  probability  of  error  for  BPSK  modulation  is  ob¬ 
tained  as 


H —  f  (1  —  $j(u))  $,,(».)  cos  (/nt)S“^  U^~  du 

X  Jo  U 


(4.12) 


provided  that  I  and  X  are  independent.  However,  this  assumption  is  discussed  in  detail 
in  Section  2.4,  where  it  is  justified  by  the  evaluation  of  the  characteristic  function  of  I  via 
two  different  methods.  Similarly,  for  OOK  modulation.  Pe  is 


,  I  f  ,1  ,  Mv  ,  ,  sin(pu)  —  sin(p  —  X)u 

+  —  /  (1  -  $/(«))  $„(») . —  - — - --du 

2tt  J0  u 


(4.13) 


where 


$,;(it)  =  exp(-~^:u2) 


S 


(4.14) 


As  shown  in  [3] -[5]  where  the  characteristic-function  method  was  applied  to  other 
systems  with  additive  interference,  expressions  similar  to  those  in  (4.12)  and  (4.13)  can 
be  evaluated  with  any  desirable  accuracy.  Therefore,  for  all  practical  purposes  results 
obtained  via  this  method  are  considered  exact  and  have  been  used  in  the  literature  for 
checking  the  accuracy  of  other  techniques.  In  particular,  the  infinite  integrals  in  (4. 12) 
and  (4.13)  are  first  truncated  (for  the  specific  acceptable  truncation  error)  and  then  the 
integration  is  carried  out  using  standard  methods  (for  the  specific  acceptable  integration 
error). 


5.  Computation  of  the  Characteristic  Function 

The  computation  of  Pe  in  (4.12)  and  (4.13)  requires  the  characteristic  function  of 
multiuser  interference  $/(u)  .  This  is  defined  as 


$/(u)  =  E[eiuI]  =  E 


iuY' 

l  — — J  m 


,(m  ) 


-sin c( - 


Ecos(  A' m  n ) 

71  =  1 


(5.1) 


where  Xmn  is  defined  in  (4.9).  We  pursue  two  approaches  for  obtaining  the  characteristic 
function. 

According  to  the  first  approach,  it  is  assumed  that  the  carrier  frequencies  of  the 
users  are  uniformly  distributed  in  a  bandwidth  of  IF.  The  phase  signature  sequences  are 
uniformly  distributed  in  the  set  of  equally  spaced  phases,  (0,2t t/M.  . . .  ,(M  —  1)2; r/M}, 
where  M  is  the  number  of  points  in  this  set.  This  means  that  the  phase  signature  sequence 
is  modulated  by  an  M-ary  phase  shift  keying  scheme.  M=2  corresponds  to  the  commonly 
used  BPSK  case. 

According  to  the  second  approach,  it  is  assumed  that  the  phase  signature  sequences 
are  continuous  and  imi/or7nh/-distributed  in  [0,2tt].  This  approximates  the  case  in  which 
the  number  of  levels  M  is  large.  The  choice  of  carrier  frequencies  distributions  is  arbitrary. 

As  it  will  become  clear  in  the  following,  the  two  approaches  not  only  correspond  to  two 
different  sets  of  assumptions  that  can  be  useful  under  a  variety  of  system  conditions,  but 
are  also  necessary  to  validate  the  use  of  characteristic-function  method  and  the  evaluation 
of  the  characteristic  function  of  other-user  interference. 
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First  we  prove  the  following  two  lemmas,  follows. 

•  Lemma  1 

Let  X  be  a  random  variable  uniformly  distributed  in  the  bandwidth  [O.TT'j.  Let 
Y  =<  X  >.  Then  the  distribution  of  Y  approaches  to  a  uniform  distribution  in 
[0,  27 r],  as  W  — co. 

•  Proof:  See  Appendix  B. 

•  Lemma  2 

Let  the  sequence  { (f>n}n  be  i.i.d.  and  uniform  in  [0,  27t]  and  let  the  sequence  {A„}„ 
be  arbitrary.  Assume  that  these  two  sequences  are  independent  for  all  n.  Then  the 
sequence  {Xn}n  defined  as 

Xn  =<  6n  +  A„  > 

is  i.i.d.  and  uniform  in  [0,  27t] . 

•  Proof:  See  Appendix  C. 


1.  Assumption  of  Uniform  Carriers 

Let  Xmn  in  (4.9)  be 

A  m  n  ~  Y  T  ra  n  T  fim  n  Y  (5.2) 

where 

7mn  =  4mn  4  kn  "F  @mn  (5-3) 

Pmn  =  u'm(n  -  1/2 )TC.  (5.4) 

{u'm}  is  assumed  to  be  i.i.d.  and  uniformly  distributed  in  a  bandwidth  IF  as  big  as  10 
THz.  Therefore,  {fimn}  for  fixed  n  are  i.i.d.  and  uniformly  distributed  in  a  bandwidth  of 
(n  —  1/2 )W/(RN),  where  R  is  the  data  rate.  For  typical  values  of  R  and  N,  this  value 
is  still  very  large  and  according  to  Lemma  1,  it  is  reasonable  to  consider  {<  ,6mn  >}  as 
i.i.d.  random  variables  uniformly  distributed  in  [0,  2tt]  .  Consequently,  Lemma  2  asserts 
that  {A'mn}  with  respect  to  m  are  i.i.d.  and  uniformly  distributed  in  [0,2tt].  Therefore, 
(5.1)  is  expressed  as 


*/(«)= 

m 


_o _ 

N 


sinc( 


cos(Xmn) 


(5.5) 
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This  is  independent  of  the  index  m  .  Hence 


$l(u)  =  {  E 


s(-V„) 


:k- d 


(5.6) 


Conditioned  on  to  in  (5.6),  it  is  easy  to  show  that  the  sequence  {cos  AT}  has  zero-mean 
and  zero  correlation.  Therefore,  this  sequence  is  a  “ p  —  mixing"  one  with  p(n)  =  0  (See  [7] 
and  [8]).  We  want  to  apply  the  CLT  to  the  sequence  {cosAra}.  To  achieve  that  we  report 
a  theorem  from  [8] . 

•  Theorem 

Let  {Aln}  be  a  second-order  stationary,  centered,  and  p-mixing  sequence;  let  a~  — >  oo 
and  £,-p( 2!)  <  oo.  Then  {sn/crn}  satisfies  the  CLT,  where  sn  —  )C”=1  A";  and  cC  = 
Vctr(sn). 

The  variance  of  the  sequence  {cos A',,}  turns  out  to  be 


2  _  j  \  +  |  cos[iu'(2n  -  l)Tc]E{cos('2dn )},  for  M  —  2 
l  2  i 


for  M  >  2 


o./ 


Since  2 9n  is  a  zero-mean  Gaussian  random  variable  of  variance  16? rjSnTc,  it  is  easy  to  show 
that  the  expectation  in  (5.7)  is  exp(— 8?r,i 3nTc).  Hence, 


2  _  (  b  +  |  cos[c</(2n  —  1)TC]  exp(  —  Stc ilnTc),  for  M  —  2 

® n  ]  { 

l  2  > 


for  M  >  2 


(5.8) 


Apparently,  all  conditions  of  the  theorem  are  satisfied  for  M  >  2.  But,  for  M  =  2.  the 


N 


sequence  is  not  second-order  stationary.  The  variance  of  the  term  ^  cos  (AT 


is 


n—  1 


2  _  /  \  +  cos[a/(2n  —  1)TC]  exp(  —  8^6nTc),  for  M  =  2 

r,  '  for  M  >  2 


(5.9) 


It  is  easy  to  show  that 


1  r 1 

lim  —  E^=1  cos[u/(2n  —  1)TC]  exp(— 8tt /3nTc)  =  /  cos(2u/T)exp(—  Svr3Tt)dt 
N  — CO  N  n 
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which  is  easy  to  compute.  Therefore,  as  N  — »  oo,  (5.9)  converges  to 


2  +  2(a2+ 62)  [(&sin(^)  -  acos(6))exp(-a)  +  a]  ,  for  M  =  2 
i  for  M  >  2 


(5.10) 


where  a  —  8i t(5T  and  b  =  2  u>'mT.  This  limit  is  positive  (non  zero)  for  almost  all  values  of  a/. 
Despite  the  fact  that  for  M  =  2  the  condition  of  stationarity  required  by  the  theorem  is  not 


satisfied,  we  have  shown  that  for  almost  any  given  uj'  the  sequence  {  jt  W  }  converges 

«,v 


to  a  positive  limit.  This  also  leads  to  the  satisfaction  of  the  condition  Yln=i  an  ~ *  °°-  as 
N  — >  oo,  stated  in  the  theorem.  We  may  thus  proceed  under  the  assumption  that  the  CLT 

N 

holds.  Then,  for  large  N,  the  term  -d=  ^  cos(.Y„)  in  (5.6)  can  be  replaced  by  (  ,  where 


n—  1 


C  is  a  Gaussian  random  variable  of  zero-mean  and  variance  a2 .  Conditioned  on  bo  and  u/, 
(5.6)  becomes 


*r(u)  =  {Ebo^  {E  [elu^stn 


bn  ‘  /  ^  T n  \  4. 

I  —§=sinc{  — 35—  K 


(K- 1) 


(5.11) 


The  inner  expectation  in  (5.11)  is  the  characteristic  function  of  a  Gaussian  random  variable. 
Therefore, 

(A'-l) 


$/(«)  -  E 


'“2  ^w{sincl-'£~TS''>)  v2 


(5.12) 


where  the  expectation  is  with  respect  to  bo  and  a/.  For  BPSK,  boe{  —  1, 1}  and  consequently 
we  obtain 


$/(u)  =  <  exp 


(  ■  ,  u'Tc  ,  ,  , 

—  (5mc(— — ))  <r- 


(  /v  —  X ) 


(5.13) 


For  OOK,  6oe{0, 1}  and  thus  we  obtain 


$l(u)  = 


2  +  2eXp 


u- 


2  N 


LO'Tr 


sinc( ) 


a2 


( A'— 1 ) 


(5.14) 


Notice  that  this  characteristic  function  does  not  depend  on  the  phase  noise  for  M  >  2. 
Thus,  when  conditioned  on  the  phase  noise  random  process,  the  pdf  of  the  other-user 
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interference  term  I  defined  by  (4.11)  and  (4.8)  turns  out  to  be  independent  of  X.  This 
provides  partial  justification  for  the  independence  of  I  and  X  in  (4.10)  and  thus  for  the 
validity  of  (4.12)  and  (4.13)  under  this  approach.  For  M  =  2,  a1  in  (5.10)  depends  on  this 
phase  noise.  However,  the  numerical  calculation  of  the  characteristic  functions  for  both 
cases  (M  =  2  and  M  >  2)  shows  a  negligible  difference  between  these  two  (see  Figure  2.6). 
Therefore,  for  M=2,  the  independence  assumption  is  also  well  made  and  consequently, 
the  validity  of  (4.12)  and  (4.13)  is  justified.  The  issue  of  the  independence  of  I  and  A' 
[necessary  for  the  validity  of  (4.12)-(4.13)]  is  raised  by  the  fact  that,  although  and  X  are 
mutually  independent  for  all  m  (and  so  are  Im  and  Imt  for  m  ^  m'),  X  and  ( Ij ,  1 2 .  •  •  •  ,  Ijy ) 
may  not  be  mutually  independent.  Final  justification  for  the  validity  of  this  approach  is 
provided  by  its  close  agreement  (refer  to  the  section  on  numerical  results,  in  particular 
to  Figures  2.9  and  2.10)  with  the  results  obtained  via  the  second  approach  that  follows. 
Under  the  second  approach,  I  and  X  are  guaranteed  to  be  independent. 


2.  Assumption  of  Uniform  Phase  Signature  Sequence 

For  the  case  in  which  the  signature  sequence  phases  are  uniformly  distributed  in  a  set 
of  equally  spaced  discrete  levels,  if  the  number  of  levels  in  this  set  is  reasonably  large,  this 
discrete  uniform  distribution  is  approximated  with  a  continuous  uniform  phase  in  [0,  2rr]  . 
Let  us  express  the  characteristic  function  of  the  interference  I  as 


*/(«)  =  Ez?  E-oEY[eluI }  (5.15) 

where  I  is  given  by  (4.11)  and  (4.8).  E y  is  the  expectation  with  respect  to  the  N(K  —  1) 
dimensional  vector  in  { Amn}  .  E-^  is  the  expectation  with  respect  to  the  A'  —  1  dimensional 

vector  in  {£>5^}.  Finally,  Es—r  is  the  expectation  with  respect  to  the  A  —  1  dimensional 
vector  u)'  =  (c o[, . . .  ,io'K).  Since  {<f>mn}  in  (4.9)  are  i.i.d.  and  uniformly  distributed  in 
[0,  27t]  for  all  m  and  n  ,  {Xmn}  are  also  i.i.d.  and  uniform  in  [0,  2tt],  for  all  m  and  n  (based 
on  Lemma  2).  Therefore, 


$l(u) 


n 


Eu 


■  u  i  (  rn  ) 

2-^-6^  -sine 

e 


COS(  Xmn  ) 


(5.16) 
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We  use  the  identity 


E 


Y 

-xmn 


JuCOs(Xmn) 


9  /'7r/2 


7T 


_  i_  f 

“  2-  J-* 

cos  ( if  cos  a; )  da*  =  Jo  ( u ) 


•  dx  — 


(5.17 


where  Jo(-)  is  the  Bessel  function  of  the  first  kind.  Using  (5.17)  the  inner  expectation  in 
(5.16)  becomes  Jo  (b^  ■  jj  ■  sine  •  which  is  independent  of  n.  Also  using  the  fact 

that  Jo(a)  =  Jq(— a)  we  obtain  the  following  expression  for  BPSK: 


(5. IS) 


Moreover,  if  we  assume  that  the  sequence  {uim}  is  i.i.ch,  (5.18)  is  independent  of  the  index 
m  and  (5.16)  becomes 


$/(«)  = 


(5.19) 


where  the  overline  is  the  expectation  with  respect  to  the  generic  to'  scattered  in  the  band¬ 
width  of  the  channel.  Similarly  ,  for  OOK  modulation.  $/(u)  is  obtained  as 


$i(u)  = 


(5.20) 


Under  the  conditions  of  this  approach,  (Ji,  J2,  •  •  ■ ,  Ik)  and  X  are  mutually  independent 
and  thus  X  and  I  are  mutually  independent;  this  justifies  the  use  of  (4.12)  and  (4.13). 


6.  PDF  of  X 

The  pdf  of  X  as  defined  in  (3.2)  depends  on  /3T,  where  2/3  is  the  Lorentzian  bandwidth 
of  the  Brownian  motion  process  A 6k(t)  . 

a)  Monte  Carlo 
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In  order  to  obtain  this  pdf  through  Monte  Carlo  simulation,  we  write  (3.2)  as 


i  F  rnS 

A  =  —  V'  /  cos  (A  dk{t))dt  (6.1) 

1  n—\  — 


where  F  is  the  number  of  divisions  of  T  into  smaller  portions  of  size  b,  where  b  —  T/F  . 
For  F  large  enough,  let 


A 9k(t)  —  A9k(n6 )  =  0„  for  (n  —  1)8  <  t  <  nS.  (6-2) 

By  substitution  of  (6.2),  (6.1)  becomes 


-A  =  cos(Q„)  =  Cr'(0)  (6.3) 

n—1 

where  {0„}^=1  is  a  sequence  of  zero-mean  Gaussian  random  variables  of  variance  Aitfinb 
and  correlation  4ir(38  min(m,  n)  .  The  statistics  of  A'  are  estimated  through  the  generation 
of  the  random  sequence  {0n}^=1  in  a  computer  by  using  (6.3).  The  probability  that  A" 
lies  in  a  set  A  —  [01,0,2],  where  ai,O2e[0, 1]  ,  is  obtained  as 

Pa  =  [  1  A{G(9)).Pe(0)cW  (6.4) 

Jrf 

where  P©  is  the  joint  distribution  of  0  =  {©i,  ©2, ©f}  and  ln[-]  the  indicator  function. 
The  Monte  Carlo  estimate  of  this  probability  is 


M 


1—  1 


(6.5) 


where  M  is  the  number  of  independent  simulation  runs  in  the  computer  and  0,-  is  generated 
according  to  the  pdf  P©.  The  sequence  {0„}^=1  generated  in  (6.5)  is  obtained  via  the 
iteration 


r  ©  1  =  1/1 

I  @ra  -  0n_i  +vn  2  <  n  <  F 


(6.6) 
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where  {vn}  is  a  sequence  of  i.i.d.  zero-mean  Gaussian  random  variables  of  the  variance 
47 r(3S.  ,It  is  easy  to  show  that  the  sequence  {©„}  has  the  required  properties. 

The  estimate  in  (6.5)  is  an  unbiased  estimate  whose  variance  is  given  by 

var(P A)  =  Pa  ajPa--  (6.7) 

The  accuracy  of  this  estimate  is  limited  by  the  statistical  fluctuations  in  the  simulation. 
More  specifically  the  accuracy  of  the  estimate  probability  of  a  quantile  q  obtained  after  M 
independent  simulation  trials  is  approximately  yj  1  /M q  .  At  low  probabilities,  the  number 
of  simulation  trials  increases  drastically,  which  prohibits  the  utilization  of  the  simple  Monte 
Carlo  method.  To  circumvent  this  problem  a  modified  Monte  Carlo  simulation,  namely 
the  Importance  Sampling  method,  is  implemented. 

b)  Importance  Sampling 

The  idea  is  that  the  input  distribution  is  biased  so  that  more  samples  lie  in  the  region 
of  interest  at  the  output.  Let  us  write  (6.4)  as 


Pa  =  [  i A[G(d)].w(0).P^(d)dd  (6.S) 

Jrf 

where  W{(P)  =  p^;^  and  P©«  is  the  biased  distribution  of  0.  The  Importance  Sampling 
estimate  of  Pa  in  (6.8)  is  given  by 


1  M 

px  = 


where  0*  is  generated  from  the  distribution  Pq * .  The  variance  of  this  estimate  is  given 
by  (see  [10]) 

W  -  P\ 

var(P*A)  =  . —  -4  (6.10) 

where  the  average  weight  W  is 

W=  f  lA[Cr(0)}W(e)Po(0)d9.  (6.11) 

Jrf 
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The  fundamental  issue  in  Importance  Sampling  is  the  determination  of  the  biasing  dis¬ 
tribution  P©*,  so  that  the  variance  of  the  Importance  Sampling  estimator  is  minimized. 
The  proposed  biasing  distribution  is  obtained  by  shifting  the  original  input  samples  Q„  by 
±cvn<5  ,  i.e., 

0^  =  0,i  +  ZanS  (6.12) 


where  Z  takes  values  in  {  —  1, 1}  with  probability  1/2  .  The  parameter  a  >  0  is  the  identity 
which  should  be  optimized  accordingly.  Given  P©  and  using  (6.12)  we  can  easily  show 


that 


W{9)  = 


Pq(£) 
Per  (9) 


exp( 


(aT)2  x 
8tt i3T  > 


cosh( 


aT  Of  \ 
4tt3T  ) 


(6.13) 


This  is  the  same  expression  which  was  given  in  [11]  in  a  similar  context.  However,  it  was 
not  proposed  in  [11]  how  to  optimize  (6.13)  in  terms  of  a. 

In  order  to  reduce  the  variance  in  (6.10)  with  respect  to  the  one  in  (6.7)  we  need  to  make 
sure  that  W  <  Pa-  In  what  follows,  although  we  do  not  provide  the  optimum  value  for 
a  such  that  this  condition  holds,  but  we  argue  for  an  approximation  in  the  hope  that  the 
numerical  results  will  provide  encouraging  results. 

First  let  rewrite  expressions  (6.4)  and  (6.11)  as 


Pa  =  !  Pe(0)cW 
J  B 


(6.14) 


PF 


=  [  W(8).Pq(9)cW 
Jb 


(6.15) 


respectively,  where  B  is  the  inverse  image  of  the  set  A,  i.e.,  B  =  G'_1(.4).  Also  let  the 
set  A  be  A  =  {a}, ae[— 1, 1],  which  is  a  reasonable  assumption  due  to  the  probabilities  of 
bins  of  small  sizes.  P@_(6)  is  the  joint  pdf  of  a  multivariate  Gaussian  with  zero  mean  and 
covariance  S,  where 


/I 

1 

1  .. 

1 

2 

2 

2 

S  =  4  tv/38 

1 

2 

3  .. 

.  3 

\1 

2 

3  .. 

•  FJ 
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It  is  easy  to  show  that 

Pq(8)  =  (27r)^.|S|^1.e_2,r/J^i+(e2_0l)2+-+<5F_^-l)21.  (6.17) 

At  this  point,  we  attempt  to  minimize  W(8)  in  (6.15)  with  respect  to  a  for  the  values 
of  8_  in  B  which  maximize  P&(8).  In  other  words,  the  effect  of  the  integrand  in  (6.15)  is 
reduced  when  it  is  maximum.  As  we  proceed,  we  find  the  optimum  value  of  a  in  (6.13), 
which  minimizes  W(8)  when  Pq_(8)  is  maximum,  and  we  also  express  the  conditions  under 
which  this  minimum  value  is  less  than  one. 

First,  (6.17)  is  maximum  when  8  is  in  the  set  S  —  {#:  8\  —82  =  . . .  =  dp}  C  PF ■  This 
reduces  (6.17)  into 

P0(^)  =  (27r)^.|S|^.e-2^,  for  8eS.  (6. IS) 

Second,  we  need  to  maximize  (6.18)  with  the  constraint  8eB.  In  other  words,  we  should 
maximize  (6.18)  for  8eS  fl  B.  But  it  is  obvious  that 

S  H  B  =  {8\  81  —  82  =  • . .  =  8p  =  2A'tt  ±  cos-1  (a),  for  Iv  integer}  (6.19) 


and  the  values  which  maximize  (6.18),  and  consequently  (6.17),  are  in  the  set 


L  =  {6:61  =8,  =  ...  =  8f  =  ±cos_1(a)}  c(SnB)cB.  (6.20) 


Next,  we  minimize  W(8)  in  (6.13)  with  the  constraint  8eL.  i.e., 


min {W{8)  - 


exp( 


iocT)-  , 
8t r/3T  > 


cosh( 


ctTOp  \ 
4t t0T  > 


8eL } 


(6.21) 


or 


mm{W{6)  = 


eMp^r) 

cosh(iSr) 


(6.22) 


where  <j>  —  cos  1(a).  A  simple  derivative  of  W(8_)  in  (6.22)  with  respect  to  a  renders  the 
necessary  condition  for  the  minimum  points.  This  is 


V  =  tanh(7y) 


(6.23) 
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where  y  =  ^  and  7  =  4^7-  For  7  <  1,  (6.23)  lias  only  one  solution  at  y  =  0  .  This  means 
that  W (9)  is  minimized  at  a  =  0  and  W{9_)  =  1  .  This  is  the  useless  trivial  reduction 
of  the  problem  to  the  simple  case  of  Monte  Carlo.  For  7  >  1,  (6.23)  has  three  roots.  At 
a  =  0  W(9)  is  maximized;  at  the  other  two  symmetric  roots  (±amjn),  W(8)  is  minimized 
and  this  minimum  value  is  less  than  one.  In  other  words,  the  minimum  of  W{0}  in  (6.22) 

for  nonzero  a.  is  achieved  when  ,  ^,)rr  >  1  and  this  minimum  is  less  than  one.  However,  for 

■iTT  pi 

practical  values  of  <j>  and  /3T,  this  condition  is  always  satisfied. 


7.  Asynchronous  System 

In  this  section,  the  analysis  developed  for  the  synchronous  system  is  extended  to  the 
asynchronous  case.  The  time  delay  rm  in  (2.1)  is  considered  to  be  a  uniformly  distributed 
random  variable  in  [0,  T],  which  represents  the  m-th  user’s  time  delay.  At  the  A’-th  receiver, 
the  matched  filter  is  synchronized  with  the  k- th  signal,  i.e.  77  =  0.  Therefore.  (4.10)  is 
still  valid  and  the  average  probability  of  error  is  given  in  (4.12)  and  (4.13)  for  BPSIv 
and  00K  modulation,  respectively.  Evaluation  of  the  characteristic  function  of  multiuser 
interference  for  the  asynchronous  case  follows  next.  Let  us  write  Tm  as 

Tm  =  (-m.Tc  4-  r'm  (7.1) 


where  tm  =  is  a  discrete  random  variable  taking  values  in  {0,1....,  Ar  —  1}  with 

equal  probability.  T'm  is  uniformly  distributed  in  [0,TC].  The  counterpart  of  (4.3)  becomes 


bm(t  —  T  ^ei[UJ'mt~u’'mrm+4>m(t-TTn)-4>k(t)  +  AOrn{t)j 

The  integral  in  (7.2)  can  be  written  as 


(7.2) 


(7.3) 


where 


rT m  r(n— IIT+tL  rnTc 

JO  ^  J (n—l)Tc  +  5i  i"-D7+4 


(7.4) 
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~T  N  .{n-1)Tc  +  r'm  rnTc 

rm  n-em+ 2  “'("“DT’c  >i=en,  +  i  ("~1)Tc  +  r™ 


By  using  (7.3),  (7.4)  and  (7.5)  in  (7.2)  irn  becomes 


i-(n-l)Tc+r^ 


'(n-l)Tc 


ei9'»"lt)dt 


<  (n-\)Tc  +  T'm 


where 


9mn(t)  —  ^  ^ mT™  4“  ^mfn  — 1)  fikn  4"  @ mn 

9mn(b)  ^ m b  Bn  d~  <Pnin  @kn  d"  @rrvn 

emn  =  b-l  •  l[14m  +  l](«)  d-  to'0  •  1[£ m+2jV](«) 

«m»  =  W  '  1[l.<m](n)  +  C}  •  l[fm  +  l,,V](«) 

.  ,  ,  A  r 1  I  <n  <  j 

1Ei,j](«)  —  \  o  other 


The  counterpart  of  (4.8)  is  obtained  as 


Im  —  7?e{dm}  —  -T  ^  ^  [emnam  ^OS  X  in  n  d"  ^mn  am  COS^inn, 


where 


+  A  Bn  •  _  { 

aZ,  —  — smc  - 

m  nr  \  0 


-A,  r'm  _Wm{Tc  -  Bn)' 
“m  =  (!  -  7fr)stnc  . — - 

C 


^"mn  =B  u'm(n  —  1)TC  H - ™--  —  —  UJmTm  +  (pm(n-l)  4>kn  d-  $mn  > 


*m»  =<  «m("  "  l/2)Tc  + 


m  -  w'  Tm  +  Omn  -  bfcn  d-  #mn  >  • 


At  this  point,  we  make  two  assumptions  similar  to  those  of  the  synchronous  case,  two 
assumptions  as  follows. 
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1.  Assumption  of  Uniform  Carriers 

As  before,  we  observe  that,  conditioned  on  r  =  (ri .  79, . . . ,  Tk+i . r/x- ) ,  { A  +  () } 

and  {A~n}  are  i.i.d.  with  respect  to  m  and  uniformly  distributed  in  [0,  2tt] .  Therefore, 


l.  m 


^  W  (  ^  m  7  '  n  —  ^  n  ^  /  n  _  ^  ^  m  n 


COS  A',' 


(7.17) 


The  sequences  {eJ^nn  cos  X+xn}n  and  {e~  n  cos  X~n}n  are  zero-mean  and  have  zero  correla¬ 
tions.  Therefore,  they  are  p-mixing  sequences.  Also  notice  that  these  two  sequences  are 
uncorrelated  with  respect  to  each  other.  By  conditioning  on  u>'m  and  by  using  the  C’LT 
(7.17)  yields 


*/(«)  =  j5vmjEw'm2r 


>7 


MS] 


where  and  r]m  are  zero-mean  Gaussian  random  variables  of  variances  a2+  and  a2_  . 


rn 


respectively.  It  is  easy  to  show  that,  for  BPSK  modulation. 


&  4-  \  ■* 

Vrn  I  1 


5  +  W  £,=,  cos(2  *'m{j  -  l!,n  -m-  M  =  2 

M  >  2 


2  ’ 


(7.10) 


a2_  =J  i  +  A/  =  2 

(?m  1  1 

2  ’ 


-1/  >  2 


(7.20) 


The  two  summations  in  (7.19)  and  (7.20),  as  N  — ►  oo.  are  equal  to  the  integral 


cos 


L  -L  L  171  I  ’ 

Omt  0m  ^  —U)mTm 


exp  (—at)  clt 


(7.21) 


where  bm  =  2uj'rnT  and  a  =  8tt0T.  Let  define  ?/m  as 


Vm  =  O-lnUm  +  amVr 


7.22) 


The  variance  of  this  zero-mean  Gaussian  random  variable  is 


Vm 


(“m)2^+  +(0im)2<Jl-  =  K)2+(«mf]^(f™) 


.23) 
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where  cr^(n)  =  \  for  M  >  2,  and  for  M  —  2, 

<?mH  =\  +  b'[a2~^  fjT  )  acos(^rb™  +u'mT'm)  +  hm  +  En  T  )  + 

77-  72  \ 

bm  sin(6m(l  -  — )  -  uj'mr !„)  -  acos(ftm(l  -  exp(-«j 


If  we  substitute  (7.22)  in  (7.18)  and  perform  the  inner  expectation,  we  obtain 

$/(«)  =  Er  |]^J  Eu 


7.24) 


e  2jV  '''•Jr 


(.2D 


r  is  an  i.i.d.  and  uniform  random  vector.  Hence,  the  outer  expectation  is  interchanged 
with  the  product  and  the  result  is  independent  of  the  index  rn . 

(A— 1) 


r 

—  2  " 

) 

yErE^' 

eTTT-^r, 

l 

i 

$l(u)  = 


The  expectation  with  respect  to  r  is  reduced  to  t'  and  L  Hence,  finally 


7.26) 


*/<»)  =  Ux>»<“ 


( A’— 1 ) 


H.27) 


n= 1 


where 


$„(u)  =  Et,E^ 


-  u“  2 
e~2N~ar, 


7. 28) 

7.29) 


0-^  =  [(cv+)2  -f-  (Of  )2]  a2{n). 

For  00K  modulation,  (7.27),  (7.28)  and  (7.29)  still  hold  and  only  a2(n)  is  divided  by  a 
factor  of  2. 


2.  Assumption  of  Uniform  Phase  Signature  Sequence 

Since  the  two  sequences  {<f>m(n- 1)}  in  (7.15)  and  {4>mn}  in  (7.16)  are  i.i.d.  and  uni¬ 
formly  distributed  in  [0,27r],  for  all  m  and  n,  the  two  sequences  {A'+n}  and  (A'~n  }  are 
also  i.i.d.  and  uniformly  distributed  in  [0,27r],  for  all  m  and  n  (Lemma  2).  Therefore. 


$/(«)= n  e 


j  T  (c 


m  EjEl  em-COS 


A'+  +«: 


El 


emn  COS-U,., 


(7.30) 
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Moreover,  the  two  sequences  {X^n  }  and  {„Y~n}  are  independent  of  each  other.  To  prove 
this  we  need  to  show  that  the  two  random  variables  X^in  and  X~  are  independent  for 
all  m,n,i  and  j.  For  all  cases,  except  the  case  in  which  i  =  m  and  j  =  n  —  1,  the  i.i.d. 
and  uniform  phase  sequences  {<f>m(n-i)}  in  (7.15)  and  {d>jj}  in  (7.16)  provide  the  proof 
in  a  straightforward  manner.  For  the  case  i  =  m  and  j  —  n  —  l,Y'+n  and  X~  are  still 
independent,  because  <f)kn  and  (j)k(n-i)  are  independent  and  uniform  for  this  case.  Using 
these  facts 


§l(u)  =  Em 

m 


N 

n 

„  n=l 


Ex+ 

m  r 


e  t 


v  + 


Ex: 


e  ? 


(7.31) 


where  the  expectation  Em  is  with  respect  to  the  all  random  variables  with  index  m.  Upon 
substitution  from  (5.17) 


$l(u)  = 

m 


+  e+ 

m  ^  m  n 


(7.32) 


For  BPSK  modulation,  e^n  and  emn  belong  to  {  —  1,1}.  Therefore,  the  inner  product  is 
independent  of  n  and  (7.32)  becomes 


(7.33) 


Moreover,  by  assuming  that  the  sequences  {u/m}  and  {r'm}  are  i.i.cb,  the  expectation  above 
becomes  independent  of  the  index  m  and,  finally,  the  characteristic  function  of  the  mul¬ 
tiuser  interference  for  BPSK  modulation  is  obtained  as 


$/(u)  = 


U  T 


Jo[N'  TrSmC 


u 


.7°  (  -  (  1  -  -  ' 


T 

1 

eT 

3 

i—.. _ 

\] 

N\ 

2  J 

)\ 

J 

(A'-l) 


7.34) 


where  the  expectation  in  (7.34)  is  with  respect  to  ad  and  r',  which  are  distributed  in 
[— TU/2,  W/2]  and  [0,TC],  respectively. 
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The  evaluation  of  the  characteristic  function  of  the  multiuser  interference  for  OOK  in 
this  case,  is  more  tedious  than  for  BPSK .  We  know  that  lm  is  a  uniform  discrete  random 


variable  which  takes  values  in  {0, 1, . . .  ,  N  —  1}.  By  assuming 


A  ,  U  _j_  . (m) 

=  Mfamb-1 


.  By  assuming 

it”') 

(7.35) 

(7.36) 

)Jo  (frf’) 

(7.37) 

(7.32)  becomes 


m  t  Ln=l  J  Lra=£m+2 
Since  cm  ,  dm  and  vm  are  independent  of  n,  then 


=  J J  Cm  [J  d 


§ I (u)  =  {  <-’>»  Cm"  d] 


S  —  /  —  I 

m 


Performing  the  expectation  above  with  respect  to  im ,  yields 


By  taking  expectation  with  respect  to  ,  b(urn>)  and  rearranging  the  terms  the  final 
result  becomes 

$/M  =  { j  +  ji TTTWTWTTb 


l[,„ . „  ,  ,  ,s, 

+  PM"  -T  )  +  *■<"  -T'  - r  )] 


(  h'  —  l) 


(7.41) 
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where 


Ru(w\t')  —  Joi^r'  sinc(—~)). 


x o'  t' 


(7.42) 


8.  Numerical  Results 

Figures  3a  and  3b  illustrate  the  pdf  for  the  random  variable  X,  which  was  defined 
in  (3.2).  According  to  Section  2.5,  by  using  the  importance  sampling  technique  this  pdf 
is  obtained  for  different  values  of  the  parameter  (3T .  As  expected  from  (3.2),  for  lower 
values  of  phase  noise  (low  f3),  this  pdf  tends  to  accumulate  around  the  point  X=1  which 
corresponds  to  the  noiseless  case.  By  using  the  importance  sampling  method  we  were  able 
to  find  the  tail  of  this  pdf  down  to  10-9. 

Figure  4  shows  the  performance  of  the  single-user  system  for  the  BPSK  case.  Degra¬ 
dation  of  this  performance  is  due  to  two  factors,  additive  white  Gaussian  noise  and  phase 
noise.  For  values  of  j3T  higher  than  0.05,  the  average  probability  of  bit  error  is  higher  than 
10-3,  even  for  values  of  Eb/No  up  to  20  dB,  where  Eb  is  the  average  bit  energy.  Reduction 
of  the  phase  noise  by  50%  (/3T  =  0.025)  reduces  Pe  to  the  range  of  10-'3.  Even  higher 
performance  is  obtained  for  the  lower  values  of  /3T.  The  lower  bound  is  achieved  when 
/ 3T  =  0  (no  phase  noise). 

In  Figure  5,  the  performance  of  the  two  schemes,  BPSK  and  OOK,  for  the  single- user- 
case  are  compared.  For  the  no-phase-noise  case,  OOK  is  worse  than  BPSK  for  3  dB.  This 
value  increases  adversely  due  to  phase  noise.  As  shown,  for  /3T  =  0.015,  the  two  curves  for 
BPSK  and  OOK  diverge  even  further  for  higher  values  of  SNR.  This  asserts  the  inferior 
performance  for  OOK  due  to  phase  noise. 

Numerical  calculation  of  the  characteristic  function,  for  typical  values  of  the  system 
parameters,  reveals  that  the  binary  and  M-ary  phase  signature  sequence  in  the  first  method 
yields  values  that  are  very  close  to  each  other.  This  not  only  validates  the  average  prob¬ 
ability  of  bit  error  formulation  for  the  case  of  binary  phase  signature  sequence  (M=2), 
but  also  suggests  considering  the  M-ary  case  (M  >  2)  (which  is  computationally  easier  to 
obtain)  as  the  representative  case  for  all  numerical  results.  Figure  6  shows  how  close  these 


two  cases  are. 


Multiuser  performance  of  the  asynchronous  system  for  BPSK  is  illustrated  in  Figure 
7.  For  5000  high  data  ra.te  users  of  data  rate  10  Mbs,  the  average  probability  of  error  is 
obtained  for  different  values  of  /3T  and  SNR.  For  this  high  number  of  high  data  rate  users 
and  for  values  of  ftT  such  as  0.015,  the  average  probability  of  bit  error  in  the  range  of 
10~6  becomes  feasible  by  introducing  more  signal  power  (higher  SNR).  The  total  optical 
bandwidth  W  and  the  spreading  gain  N  are  set  to  be  10  THz  and  1000,  respectively.  The 
multiuser  capability  of  the  same  system  is  obtained  in  Figure  S.  For  a  Pe  in  the  range 
of  10-4,  this  system  can  support  almost  10000  high  data  rate  users  by  tolerating  a  phase 
noise  as  high  as  /3T  =  0.015,  while  keeping  Eb/No  higher  than  12  dB. 

Performance  of  the  synchronous  and  asynchronous  systems  using  the  two  proposed 
methods  of  evaluation  of  the  characteristic  functions  is  obtained  for  BPSK  and  OOIv  in 
Figure  9  and  Figure  10,  respectively.  For  the  proposed  values  of  parameters  indicated,  the 
performance  of  the  synchronous  system  using  cither  method  and  the  asynchronous  system 
using  the  first  method  are  exactly  the  same.  In  fact,  we  obtained  computationally  exactly 
the  same  values  for  the  corresponding  characteristic  functions.  The  asynchronous  system 
using  second  method  is  slightly  different  from  the  others.  Also  notice  that,  in  the  OOIv 
case,  the  degradation  of  the  performance  is  mainly  due  to  the  single-user  contribution  of 
error  in  the  total  average  probability  of  error. 

The  close  agreement  of  the  results  provided  by  the  two  methods  observed  in  Fig¬ 
ures  9  and  10  justifies  the  use  of  the  first  method,  which  corresponds  to  more  realistic 
system  modeling  assumptions,  whereas  the  second  method,  which  is  analytically  very  ac¬ 
curate  (actually  exact  within  the  framework  of  the  characteristic-function  method),  relies 
on  modeling  assumptions  that  are  less  realistic.  Therefore,  depending  on  the  modeling 
assumptions  we  can  use  the  most  suitable  method  with  full  knowledge  of  the  accuracy  of 
the  method. 

9.  Conclusions 

In  this  paper  we  provided  the  modeling  and  the  tool  for  analyzing  the  performance 
of  an  RC  CDMA  network  with  coherent  optical  detection.  The  average  bit  error  prob¬ 
ability  of  this  multiple-access  scheme  was  evaluated  using  the  characteristic  function  of 
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multiuser  interference  at  the  output  of  the  matched  optical  filter.  Based  on  two  set  of 
assumptions,  we  proposed  two  approaches  for  obtaining  the  characteristic  function.  These 
two  methods  were  numerically  evaluated  and  the  results  proved  to  be  very  close  to  each 
other.  Both  phase  noise  and  thermal  noise  were  taken  into  account  in  the  analysis.  The 
effects  of  AWGN  and  phase  noise  were  studied  numerically  for  single-user  systems  with 
spreading  and  despreading  and  for  multiuser  systems.  Importance-sampling  techniques 
were  developed  to  simulate  a  given  function  of  phase  noise.  It  was  shown  that  thousands 
of  high  data  rate  users  can  communicate  reliably  in  a  multiaccess  optical  channel  corrupted 
by  AWGN  and  phase  noise.  Time-synchronous  as  well  as  asynchronous  systems  were  an¬ 
alyzed  in  this  context.  The  nearly  equal  performance  of  synchronous  and  asynchronous 
schemes  indicated  that  synchronization  does  not  enhance  performance  significantly.  Bi¬ 
nary  phase-shift-keying  (BPSK)  and  on-off-keying  (00K)  data  modulation  schemes  were 
considered.  Considerable  performance  degradation  of  00K  due  to  phase  noise  left  the 
BPSK  scheme  as  the  dominant  viable  scheme. 

The  performance  evaluation  of  RC  CDMA  established  the  potential  advantage  in 
employing  hybrids  of  wavelength-division  multiple-access  (WDM A)  and  CDMA  to  combat 
inter-carrier  interference  in  dense  WDMA  systems. 
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Appendix  A 

In  this  appendix,  we  provide  the  average  probability  of  bit  error  for  single  and  mul¬ 
tiuser  systems.  BPSK  and  00K  modulation  are  both  considered. 

1.  Single  User  Analysis 

The  probability  of  error  Pe  for  BPSK  modulation  is 


Y  >  p\/P\b[k)  =  -1 


+  -^Pr 


Y  <  py/P\b[k)  =  1 


(A.  1) 


where  py/P  is  the  threshold.  Upon  substitution  for  Y  from  (3.1) 


pe  =  \Pr  [v  >  p  +  X]  +  ]-Pr  [;/  <  p-  X] . 


(A.2) 


The  random  variable  X  takes  values  in  [—1,1].  Pe  in  (A. 2)  takes  two  possible  forms, 
depending  on  the  values  of  A', 


pe  =  ,  i  -  Yr  [p  -  X  <  ,  <  p  +  X/X  >  0]  w.p.  p>  =  Pr  [X  >  0]  3) 

\  +  \Pr  [p  +  A'  <  T]  <  p  —  X/X  <  0]  w.p.  q*  =  Pr  [X  <  0] 


where  “w.p.”  means  “with  probability”.  By  taking  the  average  of  (A.3),  we  obtain 

Pe  =  ^  +  Y Pr  [p  +  X  <  Ij  <  p  -  X/X  <  0] 


p 


—Pr  [p  —  X  <  r)  <  p  +  X/X  >  0]  . 


(A. 4) 


Next  we  find  the  two  probabilities  in  (A. 4).  The  first  one  is 

—  Pr  [p  +  A'  <t ]  <  p  —  X.  X  <  0] 


< 1 *  J- 1 


Pr  [p  +  x  <  i]  <  p  —  x\  fx(x)dx 


A.  5) 


where  fx(~)  is  the  pdf  of  the  random  variable  X.  By  using  the  functions  <£(.)  and  Q{ . ) 
(A. 5)  takes  the  form 


i  r° 

q*  J- 1 


$  (p  —  x) 


2 PT  \  x  I2PT 

A  Q  I  ( p  +  x )  i 


N„ 


No 


fx(x)dx  (A. 6) 
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where  Q(-)  and  $(•)  are  related  to  the  standard  normal  distribution  as  follows: 


Q(°) 


r  > 

Ja  \j2x 


"*2/ 2  •  dx 


$(a)  =  [  — =e  1  '  dx 

J-  oo 


Similarly,  the  second  integral  in  (A. 4)  is 


P*  Jo 


1  —  $  (p  —  x) 


2PT 


No 


Q  ip  +  x) 


2  PT 


No 


fx(x)dx.  (A.  7) 


Upon  substitution  of  (A. 6)  and  (A. 7)  in  (A. 4)  the  final  result  is  obtained  as 


Pe  =  x 


(US) 


The  overlines  in  (A. 8)  indicate  expectation  with  respect  to  the  random  variable  A".  Simi¬ 
larly,  for  00K  modulation  Pe  is  obtained  as 


Pe  = 


[A.  9) 


It  is  easy  to  show  that  the  optimum  value  of  p  for  (A.  8)  and  (A.  9)  is  0  and  1/2.  respectively. 

2.  Multiuser  Analysis 

The  probability  of  error  Pe  for  BPSK  modulation  is  obtained  from  (A.l)  as  follows.  If  we 
replace  “77”  with  11 1  +  77”,  (A. 2), (A. 3),  and  (A. 4)  are  still  valid  in  the  multiuser  case.  The 
first  probability  in  (A. 4)  is 

—  Pr  [p  +  X  <  I  +  77  <  p  —  X,X  <  0] 


1 

—  Pr[p  +  x  <  I  +  77  <  p-  x]  fx(x)dx. 

9  J- 1 


A. 10) 
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Here  we  made  the  assumption  that  I  and  X  are  independent.  The  validity  of  this  as¬ 
sumption  is  justified  in  Section  2.4  by  evaluating  the  characteristic  function  of  multiuser 
interference.  Upon  substitution  of  Pr  [•]  in  (A. 10)  by  an  integral,  we  obtain 


,  ,0  j-p-x 

-  /*(«)/ 

(1  J- 1  J  p+x 


fi+ri{y)  dydx 


.4.11) 


where  fi+v(-)  is  the  pdf  of  I  +  rj.  Since  the  pdf  //+,,(•)  is  a  real  function,  it  is  easy  to  show 
that 


fl+rj(y )  =  -  /  [Re  {§i+n(u)}  •  cos (yu)  +  Im  {<j>i+v(u)}  sin(yu)]  du 

7T  /n 


1.4.12) 


where  i?e{-}  and  Im{-}  denote  the  real  and  imaginary  parts,  respectively,  of  their  complex 
arguments.  Considering  that  I  and  p  are  independent  and  that  $/(u)  turned  out  to  be 
real  (see  Section  2.4), 


//+„(!/)=-  /“$/(<.)  #, 
^  Jo 


(it)  cos (yu)du. 


[-4.13) 


Substituting  (A. 13)  in  (A. 11)  and  performing  the  integration  JJ+*,  the  first  probability  in 
(A. 4)  is 

2  r°°  ir°  1  du 


$/(u)  <h,;('ti)  cos {pu) 


fx  ( x )  sin(ua’)  dx 


*r  Jo 

Using  the  same  procedure,  the  second  probability  in  (A. 4)  is 


;  .4.14) 


*P  Jo 


$/(u)  Qv(u)  cos  (pu) 


fx  ( x )  sin(u.r)  dx 


du 

u 


(.4.15) 


Combining  (A. 14)  and  (A. 15)  in  (A. 4)  we  obtain 


1  1  ^  ^  sin(uA’)  , 

Pe  =  - /  <&/(u)  §n(u)  cos  (pu) - ; -  du , 

2  7 T  J0  U 


where 


sin  (uX)  =  J 


fx (x)  sin( u x )  dx. 


(.4.16) 
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To  evaluate  Pe  for  00K  modulation  we  follow  the  same  method  as  before  and  obtain 


11  roo  sm(  pu )  —  sin(p  —  X)u 

=  o  -IT  \  *'(«)  $„(“)  '  — — - — - —  '  du. 

2  2tt  J0  u 


.4.17) 


In  order  to  put  (A. 16)  in  a  more  meaningful  format  and  also  to  facilitate  the  computations, 
we  rewrite  (A.  16)  as 


Pe  = 


1  i  r 

2  7T  Jo 


sin(uA') 

$rj(u)  COS  (pu) - U.U 


1  f00  .  sin(u.A) 

4 —  /  (l-$/(ti))  $,,(u)  cos  (pu) — - -A/. 


7T 


.4.  IS) 


The  third  term  in  (A.  18)  includes  the  contribution  of  multiuser  interference  in  BER.  The 
first  two  terms  are  single  user  contributions.  Therefore,  by  using  (A. 9)  in  (A. IS)  we  obtain 


,2PT 


-An 


1  .  ,  .  . .  ,  sin(uA’)  , 

4 —  /  (1  —  $j(u))<f>n{u)cos(pu) - 1 - -du. 

ft  Jo  u 


44.19) 


Similarly,  for  00K  modulation,  Pe  is 


l2PT 


No 


No 


1  x  ,  .  sin( pu )  —  sin( p  —  -A ) it  . 

4-  —  /  (1  -  $/(«))  $„(«) - — - - - -<*« 

2?r  Jo  u 


(-4.20) 


where 


■r  /  \  ,  -A0  2  . 

=  exp(-^;u  ). 


:.4.2i: 


Appendix  B 

In  this  appendix,  the  distribution  of  a  random  variable  which  is  the  mod.2ir  of  a 
uniform  random  variable  is  derived.  Let  A’  be  a  uniform  random  variable  distributed  in 
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the  bandwidth  [0,  W).  We  are  interested  in  knowing  the  distribution  of  the  random  variable 
1”,  which  is  defined  as 


Y  —<  X  >  .  (B.  1) 

The  random  variable  Y  is  expressed  in  terms  of  some  disjoint  intervals  as  follows: 

y  _  f  X  —  2m:  if  2m:  <  A'  <  2 (n  +  l)7r;  0  <  n  <  M  „  9 

"  1  A'  -  2Mtt  if  2Mtt  <  A  <  IF  (  ) 


where 

M  =  LW/2ttJ  . 

For  ae[0, 27r],  the  event  (Y  <  a)  is  expressed  as 
M— l 

[J  (A  —  2n7r  <  a,  2nir  <  X  <  2 (n  +  1)tt)  [J(.Y  —  2M~  <  a,  2nir  <  X  <  14')  [B. 3) 

n= 0 

The  probability  of  (B.3)  is  the  summation  of  the  probabilities  of  the  individual  events  in 
(B.3).  By  using  the  conditional  probability  and  after  some  simplifications,  Fy{a)  (the 
cumulative  distribution  function  of  Y)  becomes 

aLW/2^J+mm(«  ,  W  -  2=r[W72TrJ) 

FYW  =  - f y -  '  t-3'4) 

x4s  W  — >  oo  ,  this  distribution  approaches  to  the  uniform  one. 


Appendix  C 

Let  {An}  be 

Xn  =<  <f) n  T  A ri  > 


(CM) 


where  {f>n}  is  assumed  to  be  i.i.d.  and  uniform  in  [0,  2rr],  for  all  n,  and  {An}  is  a  sequence 
of  random  variables  with  arbitrary  distribution.  These  two  sequences  are  independent  of 
each  other  for  all  n.  Here  we  establish  two  claims: 


•  Claim  1:  {Xn}  are  uniform  for  all  n. 


•  Proof  :  It  is  a  known  fact  that 

Ah,.  | given  A„  =  <£„.  (C. 2) 

The  pdf  of  Xn  is  obtained  as 

fxjx)  =  j  fXn \\n{x\X)f\n(X) clX .  (C. 3) 

By  using  (C.2)  (C.3)  becomes 


fxJx)  =  J  U„(x)f\nMdX  =  Unix).  ( CA ) 

•  Claim  2:  {Ara}  are  independent  for  all  n. 

•  Proof  :  We  show  that  any  two  random  variables  Xn  and  X i  [n  ^  /)  are  independent. 
To  show  this,  it  suffices  to  show  that,  for  any  function  g(  * ), 

£(9(X„)9(A',)]  =  £(9(A„)]£[!/(A',)].  (C.  5) 

Given  An  and  A /  in  (C.l),  the  first  side  of  (C.5)  conditioned  on  these  values  is 

E[g{<  <Pn  +  An  >).g{<  <pi  +  A /  >)|A„,  A/].  (C  .6) 

Since  (j)n  and  4>i  are  independent,  (C.6)  becomes 

E[g{<  4>n  +  An  >)|A„]  •  E[g{<  6i  +  A;  >)|A/]  (CAT) 


or 

E[flr(A:n)|An]-f7[</(Ah)|A,]. 

From  Claim  1: 

Xn  =  Xn  | given  An 
A '/  =  Ah  | given  At. 

Therefore,  (C.8)  is 

£MA„)]  ■  £[9(A«)]. 

This  means  that 

E[g(Xn)g(Xi)\\n,  A,]  =  £[</(A„)]  ■  E[g(Xi)}. 


(CAS) 

(C'.9) 

(C'.IO) 

(CA11) 


Taking  expectation  with  respect  to  A„  and  A;  gives  (C.5)  and  this  completes  the  proof. 
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Fig.l  Coherent  detection  of  the  optical  signal  S(t) 


Fig  2  Detection  of  the  electrical  signal  T[ (t) 


Fig.  6.  Comparison  of  binary  and  M-ary  phase  signature 
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